The outer crust of non-accreting cold neutron stars 
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The properties of the outer crust of non-accreting cold neutron stars are studied by using modern 
nuclear data and theoretical mass tables updating in particular the classic work of Baym, Pethick 
and Sutherland. Experimental data from the atomic mass table from Audi, Wapstra, and Thibault 
of 2003 is used and a thorough comparison of many modern theoretical nuclear models, relativistic 
and non-relativistic ones, is performed for the first time. In addition, the influences of pairing and 
deformation are investigated. State-of-the-art theoretical nuclear mass tables are compared in order 
to check their differences concerning the neutron dripline, magic neutron numbers, the equation of 
state, and the sequence of neutron-rich nuclei up to the dripline in the outer crust of non-accreting 
cold neutron stars. 



I. INTRODUCTION 

Neutron stars can be observed as pulsars by their 
light house effect. The matter in non-accreting cold 
neutron stars is in its ground state in nuclear equilib- 
rium which means that the energy cannot be lowered 
by strong, weak, or electromagnetic interactions. Matter 
in equilibrium concerning weak interactions is termed (3- 
cquilibrated matter or matter in /^-equilibrium. Stars 
are bound by gravity and have to be charge neutral, oth- 
erwise they would be unstable and explode because of 
repulsive Coulomb forces. 

Neutron stars consist of an atmosphere of electrons, 
nuclei, and atoms. Only a fraction of the electrons are 
bound to nuclei. The ground state of the nuclei in 
this regime, with a mass density of p < 10 4 g/cm 3 , is 
56 Fe. The equation of state was calculated by Feynman, 
Metropolis, and Teller 

In the context of neutron stars with temperatures 
above typically 100 eV, a liquid layer is present between 
the atmosphere and the solid crust due to e.g. hydro- 
gen/helium burning, where nuclei and electrons are in 
a liquid phase called the "ocean" @|. However, in this 
work, we focus on non-accreting neutron stars at zero 
temperature which is a good approximation for cold neu- 
tron stars. 

One assumes complete ionization of the atoms, when 
the spacing between nuclei becomes small compared to 
the Thomas-Fermi radius ttf — o-qZ^ 1 ^ of an isolated 
neutral atom. In this equation, ao is the Bohr radius 
and Z the charge number. The mass density approxi- 
mately amounts to p ~ Am u nisr, where A is the mass 
number, m u the atomic mass unit, and the num- 
ber density of nuclei which depends on the radius of a 
spheric nucleus whose volume is the average volume per 
nucleus, 47rr 3 /3 = 1/njy [3j- By combining the last three 



equations, one finds that the outer crust of cold neutron 
stars begins when p ~ 10 4 g/cm 3 3> 3AZ g/cm 3 . This 
shell consists of nuclei and free electrons. The equation 
of state was originally calculated by Baym, Pethick, and 
Sutherland (BPS) g. The BPS model is valid for zero 
temperature (T — 0) which is a good approximation for 
the crust of non-accreting cold neutron stars. We will 
describe their model of this mass-density regime in more 
detail in the next paragraph. The inner crust of neu- 
tron stars begins when neutrons start to drip out of the 
nuclei at p ~ 4.3 x 10 11 g/cm 3 . This happens because 
the equilibrium nuclei become more and more neutron- 
rich, and finally no more neutrons can be bound to nu- 
clei. At p ~ 2 x 10 14 g/cm 3 nuclei do not exist anymore, 
signalling the end of the neutron star crust. The equa- 
tion of state of the inner crust was calculated by Baym, 
Bethe, and Pethick [5j and another equation of state of 
this regime was derived by Negele and Vautherin [||. 
Also a relativistic mean field model has been used to de- 
scribe the density regime of the neutron star crust within 
the Thomas-Fermi approximation (see and references 
therein). For higher densities, the nuclei disintegrate and 
their constituents, the protons and neutrons, become su- 
perfluid. Muons and hyperons also appear in this shell. 
The equation of state in this density regime is not well 
known, but it can be modeled by using non-relativistic 
many-body theories @ or relativistic nuclear field the- 
ories |9t llOl IllL Il2j. At extremely high densities, even 
the protons, neutrons, and hyperons disintegrate to their 
constituents, the quarks. If neutron stars with huge cen- 
tral densities exist then they can contain a quark core fl3| 
which probably is color-superconducting 
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In this paper, we focus on the outer crust of non- 
accreting cold neutron stars. It contains nuclei and free 
electrons. The latter become relativistic above p ~ 
10 7 g/cm 3 . The nuclei are arranged in a body-centred 
cubic (bec) lattice. The contribution of the lattice has a 
small effect on the equation of state but it changes the 
equilibrium nucleus to a larger mass number and lowers 
the total energy of the system because it will minimize 
the Coulomb interaction energy of the nuclei. The latter 
are stabilised against /3-decay by the filled electron sea. 
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At p ~ 10 4 g/cm 3 , 56 Fe is the true ground state. With 
increasing mass density, it is not the true ground state 
anymore because the nuclei capture electrons, emit neu- 
trinos and become neutron richer. When the mass den- 
sity p ~ 4.3 x 10 11 g/cm 3 , the so-called neutron driplinc 
is reached. Neutrons begin to drip out of the nuclei and 
become free. As soon as neutrons begin to drip out of the 
nuclei, the outer crust stops and the inner crust begins. 

The composition of the outer crust of non-accreting 
cold neutron stars was investigated by Baym, Pethick, 
and Sutherland (BPS) in a classic paper 1971 j4J. They 
calculated the equation of state and the sequence of nu- 
clei which occur in the outer crust of non-accreting cold 
neutron stars. They used the nuclear data from the 
droplet model of Myers and Swiatecki 0] . The equation 
of state of the outer crust of non-accreting cold neutron 
stars is still commonly taken nowadays from BPS al- 
though it is based on nuclear data of the mid sixties of the 
last century. Haensel, Zdunik, and Dobaczewski (HZD) 
in 1989 [13 used a Skyrme Hartree-Fock-Bogolyubov 
(HFB) calculation in spherical approximation for the pa- 
rameter set SkP 01 j hence ignoring effects from defor- 
mations, and the droplet model from Myers 0] in order 
to update the results of BPS. Haensel and Pichon (HP) 
in 1994 used the experimental nuclear data from the 
atomic mass table of 1992 from Audi and Wapstra [2l| 
and the theoretical nuclear mass tables of the droplet 
models from Moller and Nix [221 , and Aboussir et al. |23| . 
A review on the inner and outer crust of accreting as well 
as non-accreting neutron stars can be found [24j. 

In view of the previous work, it seems to be more than 
timely to reinvestigate the properties of the outer neu- 
tron star crust with up-to-date and state-of-the-art ex- 
perimental [25| and theoretical mass tables which became 
available in the last few years via the Brussels Nuclear Li- 
brary for Astrophysics Applications (BRUSLIB) |26l |2jj 
and by Dobaczewski and coworkers [2£| for Skyrme based 
models and by Geng, Toki and Meng for a relativistic 
model [29|. The nuclear models and their mass tables 
used in this work, as listed in detail in Table are taken 
to update the results of BPS, HZD, and HP. For the first 
time, the differences of various nuclear models concern- 
ing the neutron dripline, magic neutron numbers, the 
equation of state, and the sequence of nuclei in the outer 
crust of non-accreting cold neutron stars are investigated 
in detail. To our knowledge, this work is also the first 
one which uses mass tables based on modern relativis- 
tic nuclear models. Additionally, effects of pairing and 
deformation of nuclei are studied. We find that the in- 
clusion of deformations for describing nuclei is crucial in 
determining the composition of the neutron star outer 
crust. 

The results of this work rely on the (unknown) proper- 
ties of neutron-rich isotopes up to the neutron driplinc. 
The sequence of neutron-rich nuclei found in this work 
are in reach to be measured by the Facility for Antipro- 
ton and Ion Research (FAIR) [3(| at the Gesellschaft fiir 
Schwerionenforschung (GSI), Darmstadt, by TRIUMF's 



Isotope Separator and Accelerator QSAC-I I) 1 311 and by 
the Rare Isotope Accelerator project (RIA) [33 ■ Detailed 
experimental determinations of the binding energy of 
hitherto unknown nuclei towards the dripline will finally 
pin down the actual sequence of neutron-rich isotopes in 
the outer crust of neutron stars. Also, the low density 
equation of state as well as its composition serves as an 
important ingredient for low mass neutron star models 
as well as for neutron star mergers and core-collapse su- 
pernovae. 

This paper is organised as follows: In Sec. [HJ we de- 
scribe how to find the equilibrium nucleus and how to 
derive the equation of state for the outer crust of a non- 
accreting cold neutron star by using the BPS model Q. 
We use natural units (in units of MeV) and set H = c = 1 
for the equations presented in the following. We finally 
convert our results to the cgs-system in order to be able to 
compare and check the results with BPS. In Sec. IIIII wc 
describe the nuclear models and their mass tables used in 
this paper. In Sec. IIVI we present our results comparing 
the various different modern nuclear models and mass ta- 
bles used here and compare our findings to the previous 
work of BPS, HZD, and HP 0, 0, |2(j ■ We also show the 
differences of the theoretical nuclear models concerning 
the neutron dripline, magic neutron numbers, the equa- 
tion of state, and the sequence of nuclei in the outer crust 
of non-accreting cold neutron stars. In Sec.[V] we sum- 
marize our results. 



II. THE BPS MODEL 

In order to find the equilibrium nucleus and calculate 
the equation of state of the outer crust of non-accreting 
cold neutron stars by using the BPS model 3, one has 
to treat the pressure P as the independent variable. We 
start with P = 9.744x 10 18 dynes/cm 2 which corresponds 
to the mass density p ~ 1.044 x 10 4 g/cm 3 . Actually, this 
special value of the pressure is also the starting value of 
the pressure in Ref. Q| and the corresponding mass den- 
sity is in good agreement with the approximation made 
in Sec. HJ p ~ 10 4 g/cm 3 3> 3AZ g/cm 3 , the mass den- 
sity at which the outer crust begins. Because the pres- 
sure in the star is increasing continuously with decreasing 
star radius, we increase P until neutron drip is reached, 
i.e. when the chemical potential of baryons is equal to 
the neutron mass, pb = m n . With given pressure P, we 
vary the mass number A, and the charge number Z and 
solve the equation of the total pressure, 

P = P e + \w L n N , (1) 

for the electronic density n e . The highest contribution 
to Eq. is the pressure of free electrons, 

1 r ke k 4 
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where E e = (k 2 + m 2 ) 1 / 2 . The electron pressure depends 
on the electron Fermi momentum which is related to the 
electron density, 
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The lattice energy is given by 
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(4) 



It has the form of the Coulomb energy with a special 
prefactor calculated in j3^| which arises because of the 
bcc lattice. The bcc lattice constant a is related to the 
number density of nuclei by 



rijya = 2 



(5) 



The latter one depends on the number density of elec- 
trons because of the neutrality condition which has to be 
fulfilled in stars, 



Zun 



(6) 



The baryon density is related to the number density of 
nuclei, 



rib = Attn ■ 
The total energy density is given by 

E to t=n N (W N +W L )+E e 



(7) 



(8) 



The prefactor 1 /3 in Eq. originates from the fact that 
the pressure of the bcc lattice Pl is one third of the 
energy density of the bcc lattice, Pl/3 = El = Wltin- 
The energy of the nuclei is obtained by 



W N = m n (A - Z) + rripZ - bA , 



(9) 



where m n is the neutron and m p the proton mass, and b 
is the binding energy per nucleon. The energy density of 
free electrons amounts to 



E e = n e n e - P e 



where 



He 



(10) 



(ii) 



is the chemical potential of the electrons and m e the 
electron mass. The quantity to be minimized at fixed 
pressure P by varying A and Z is the baryon chemical 
potential 



gtot + P _ WN+jWL + Zjle 



lib 



A 



(12) 



This procedure has been done for different nuclear 
models. All of them contain data in tabular form for 



A, Z, or N = A — Z, and the corresponding binding en- 
ergy B or binding energy per nucleon b = B/A. We use 
the data of the nuclear models listed in Table |U 

The places at which a phase transition from one to 
another equilibrium nucleus happens can be found by 
varying the pressure as long as its difference to the pres- 
sure of the precise point of the transition becomes small. 
As the transition from one nucleus to the next one takes 
place, P and /z of both phases are equal but there will be 
jumps in fi e and n e when the proton numbers Z of both 
phases are unequal. The baryon density rib and the mass 
density 



E u 



(13) 



will jump accordingly and are approximately given by 



n b ~n b ~ n e 



A' A 



Z' ZJ ' 
Ap ^ An b ^ Z/A 
~J ~ ~ ~ Z'/A' ~ ' 

The adiabatic index T is defined by 

r= lib dP_ 
P drib 



(14a) 
(14b) 

(15) 



At the transition point, the adiabatic index jumps to zero 
because the pressure in both phases is equal. 



III. DESCRIPTION OF THE NUCLEAR 
MODELS USED IN THIS PAPER 



In this section, we describe the nuclear models used in 
this work as listed in Table [I] For comparison, we add 
the results from BPS which were derived by using mass 
tables from a droplet mass formula from the sixties. 

The most recent mass table of the finite range droplet 
model (FRDM) H3 lists 8979 nuclei ranging from 16 O 
to 136 extending from the proton dripline to the neu- 
tron driplinc. The mass table of FRDM is based on a 
macroscopic finite range droplet model including a folded 
Yukawa single particle potential. It gives so far the best 
parameterization of masses of known nuclei throughout 
the nuclear chart. 

Non-relativistic Skyrme parameterizations are the well 
known and widely used SkM* [HI, SkP [H and SLy4 
parameter sets I37J. Set SkM* originates from the pa- 
rameter set SkM [38( which is fitted to nuclear matter 
properties and properties of nuclei. Set SkM* is corrected 
for the systematically too high binding energies and too 
low fission barriers of set SkM. The set SkP derives from 
Skyrme-Hartree-Fock-Bogolyubov (Skyrme-HFB) calcu- 
lations with particular emphasis on pairing effects and 
the description of neutron-rich nuclei. The parameters 
are fitted to nuclear matter properties, in particular the 
symmetry energy and properties of 16 and 208 Pb to fix 
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the surface energy. Also set SLy4 is derived to describe 
in particular neutron-rich isotopes and the (theoretical) 
neutron matter equation of state of Wiringa et al. 
in order to improve the isospin property away from the 
/3-stability line. Nuclear matter properties, the neutron 
matter equation of state and the binding energies 
and radii of the doubly magic nuclei 40,48 Ca, 132 Sn, and 
208 Pb were utilized for the fit. Mass tables of these pa- 
rameter sets were performed by Dobaczewski et al. for 
nuclei up to Z — 108 from the proton to the neutron 
driplinc within the Skyrme-HFB approach including ef- 
fects from deformation and posted at a publicly available 
web page (see |2^j). 

The other Skyrme parameterizations are taken from 
the BRUSLIB web pages [2(| |53 and are commonly de- 
rived by fitting the masses of about 2000 known nuclei 
including effects from deformations. Different approxi- 
mations schemes have been used, however. Sets SkSC4 
and SkSC18 use the Extended Thomas-Fermi plus Struti- 
nsky Integral (ETFSI) approximation for the actual cal- 
culation of nuclei. The macroscopic part is described by 
an extended Thomas-Fermi approximation, the shell cor- 
rections included by the Strutinsky integral method, and 
the pairing energy given by the BCS approximation. The 
binding energy of about 1700 nuclei were fitted togener- 
ate mass tables up to Z = 130, 115, respectively |2f|. The 
set MSk7 originates from a Skyrme-Hartree-Fock calcu- 
lation with a ten parameter Skyrme force along with a 
four parameter <5-function pairing force using the BCS 
approximation for the pairing energy. Its mass table 
extends to Z = 120 pq. The sets BSk2, and BSk8, 
however, are generated by a full Skyrme-HFB calcula- 
tion. The binding energies of 2149 nuclei were fitted for 
these latter two sets using the experimental mass table of 
2003 25] . The corresponding mass table contains 9200 
nuclei up to Z = 120 lying between the neutron and the 
proton driplines (2?|]. The newest parameter set along 
this line, dubbed BSk9, has been constrained by fixing 
the symmetry energy to 30 MeV [2?) and is not consid- 
ered in the following. All Skyrme based mass tables used 
in this work take into account effects from deformations. 
We take the sets SLy4 and BSk8 as state-of-the-art and 
the most representatives ones for cross comparison to the 
other approaches used in this work (FRDM and relativis- 
tic models). 

Relativistic nuclear field theories used here are based 
on the exchange of mesons or relativistic point-couplings 
between nucleon fields in the mean-field and the no-sea 
approximation. The effective Lagrangians behind the 
parameter sets NL3, NL-Z2, and TMA contain the ex- 
change of scalar, vector, and isovector mesons. Sets NL3 
and NL-Z2 include scalar selfinteraction terms, set TMA 
in addition vector selfinteraction terms in the effective 
Lagrangian. The parameter set NL3 was developed in 
particular to describe isospin effects. Besides binding en- 
ergy and charge radii of 10 nuclei, neutron radii were 
included into the fit procedure . A selfconsistent mi- 
croscopic correction to the spurious center-of-mass mo- 



Model 


Set 


Comments 


Refs. 


Droplet 


BPS 


used by BPS [4] 


[16] 




FRDM 


Finite Range Droplet Model 


[34J 


Experiment 




Atomic mass table 2003 


f25] 


Non- 


BSk2 


Skyrme HFB 


[27] 


relativistic 


BSk8 


Skyrme HFB 


[27] 




MSk7 


Skyrme HF + BCS 


[26J 




SkM* 


Skyrme HFB 


[28, 35] 




SkP 


Skyrme HFB 


[28, 36] 




SkSC4 


ETFSI method + BCS 


[26J~ 




SkSC18 


ETFSI method + BCS 


[26J 




SLy4 


Skyrme HFB 


[28, 37] 




SLy4HO Skyrme HFB 


[28, 37] 


Relativistic 


Chiral 


Chiral effective model 


[40]-[43j 




NL3 


Nuclear field theory 


[44, 45] 




NL-Z2 


Nuclear field theory 


[46] 




PCF1 


Point coupling model 


[47] 




TMA 


Nuclear field theory 


[22] 



TABLE I: The nuclear models used in this paper. 



tion is performed for the set NL-Z2. Its fit encompasses 
binding energies, diffraction radii, surface thicknesses, 
charge radii and spin-orbit splittings from a total of 17 
nuclei |^. The recent par ameter set TMA [2jj updates 
the sets TM1 and TM2 which were fitted to binding 
energies and charge radii for low mass numbers (TM2) 
and high mass numbers (TM1). The fit parameters of 
the set TMA are chosen to be mass number dependent 
so as to have a good description of the properties of light 
and heavy nuclei [2!j. Note that the vector field selfin- 
teractions result in a (vector) selfenergy of the nucleon, 
which is similar to Relativistic Brueckner-Hartree-Fock 
calculations, therefore mimicking many-body effects of 
higher order beyond the usual relativistic mean-field de- 
scription (see e.g. and references therein) . The rela- 
tivistic point coupling model, developed in [43, consists 
of four, six, and eight-fermion point coupling terms in 
the effective Lagrangian. Its parameter set PCF1 is de- 
termined by a list of similar observables as in the fit for 
the set NL-Z2 (see above). Pairing effects are usually 
included by a standard <5-force within the BCS approxi- 
mation. The mass tables for the sets NL-Z2 and PCF1 
range from Z = 26 to Z = 140, the one for the set NL3 
(in spherical approximation) from Z — 20 to Z = 130 up 
to several nuclei behind the neutron dripline to look for 
particle stable neutron-rich islands of stability. Note that 
only even-even nuclei are computed for our purposes as 
only those nuclei can possibly appear in the outer crust 
in nuclear /3-equilibrium. 

The chiral effective Lagrangian used for the set de- 
noted as 'Chiral' is build on the nonlinear realization of 
the chiral SU(3)xSU(3) symmetry [U IU |H |H as mo- 
tivated from quantum chromodynamics (QCD). The full 
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nonet of scalar, pseudoscalar, vector mesons is taken into 
account as well as the baryon octet in the mean-field and 
no-sea approximation and a dilaton field. The hadron 
masses, meson masses and baryon masses, are not addi- 
tional input parameters but are generated by the vacuum 
expectation values of the scalar fields of the effective La- 
grangian (spontaneous chiral symmetry breaking besides 
explicit symmetry breaking) as dictated by the properties 
of QCD. Pairing effects are included by the BCS scheme. 
The SU(3) chiral model describes nuclear matter as well 
as properties of nuclei |4jJ, 03 > hypernuclei |4l| and neu- 
tron stars pHIT^.l43| . The mass table used extends from 
Z = 1 6 to Z = 100. 

For the relativistic models, spherical calculations with 
and without pairing effects have been performed for the 
sets NL3, NL-Z2, PCF1 and Chiral. Mass tables of cal- 
culations including effects from deformations are publicly 
available for NL3 0] and TMA [2<| and are taken for 
comparison to FRDM and the Skyrme-based parameter- 
izations. The mass table for NL3 of 0] lists 1315 even- 
even nuclei up to Z — 98 using the BCS pairing scheme 
with constant pairing gaps, while the one for TMA [2jj 
contains nuclei from the neutron to the proton driplinc 
up to Z — 100 with BCS type 5-force pairing. 



IV. RESULTS 



nuclear models, we update the results of 0, IrTL |20| . 

Besides, we study differences between the theoretical 
nuclear models directly. The location of the neutron 
dripline is of great importance for our investigations, be- 
cause their position in the nuclide chart is decisive if a 
neutron rich nucleus has a chance to be present in the 
outer crust of a non- accreting cold neutron star or not. 
If the nucleus is behind the neutron dripline, it is un- 
stable and will emit neutrons even for /3-equilibrium and 
large electron fractions. Therefore, the nucleus can not 
be present in the sequence of nuclei in the outer crust. 
Of course, this is not the case if the nucleus is located 
before the neutron driplinc. 

In this section, we subdivide our results in several cat- 
egories in order to have a better structure of the results 
in our paper. In Sec. II V Al we show the equations of 
state by using various nuclear models. In Sec. IIVBI 
we discuss the effect of pairing on the neutron dripline 
and on the sequence of nuclei. In Sec. IIV CI the non- 
relativistic parameterizations of BRUSLIB |2g, |2jj and 
Dobaczewski |28j are compared. In Sees. IIV Dl and ITY El 
we compare relativistic and non-relativistic models for 
the neutron dripline and sequences of nuclei to each 
other. In Sec. IIV Fl we show the results for the state-of- 
the-art nuclear models and discuss their proton numbers 
in the sequences of nuclei. Finally, in Sec. IIV Cfl we make 
a comparison to previous works. 



In this section, we present our results obtained within 
the model of Baym, Pethick, and Sutherland (BPS) Q 
as the equation of state and the sequences of nuclei of 
the outer crust of non-accreting cold neutron stars. Var- 
ious nuclear models as listed in Tableware used for this 
purpose. As shown in Sec. ^ the binding energy B or 
rather the binding energy per nucleon b together with 
the mass number A and the proton number Z are the 
input parameters for the BPS model, see Eq. ©. Dif- 
ferent nuclear models, of course, usually have different 
binding energies per nucleon for the same nuclei resulting 
in different sequences of neutron-rich nuclei in the outer 
crust of non-accreting cold neutron stars. Because most 
of the binding energies of nuclei with low mass numbers 
are known precisely from experiments, we prefer using 
the nuclear data from the atomic mass table 2003 from 
Audi, Wapstra, and Thibault [25| whenever possible. In 
this context, we mention that we do not take any esti- 
mated (non-experimental) data of the atomic mass ta- 
ble (2{| into account. If the corresponding nuclei are not 
listed in the atomic mass table [25], we use the data of 
theoretical nuclear models for calculating the sequences 
of nuclei which are present in the crust of non-accreting 
cold neutron stars. There is only one exceptional case: 
The nuclear data of the droplet model from Myers and 
Swiatecki 0] which is used for calculating the original 
sequence of nuclei obtained by BPS is not modified by 
nuclear data from the atomic mass table [2f| because we 
want to compare our new results with the original ones 
from BPS ||. By using the newest and most modern 



A. Equations of State 

Fig.^shows the equations of state, i.e. the pressure as 
a function of the mass density. For low mass densities, 
the nuclei appearing in the outer crust of non-accreting 
cold neutron stars are the same for all nuclear models 
because they are given by the experimental data of the 
atomic mass table j2^- Hence, the equations of state 
in Fig. n are almost the same. Only the set BPS ex- 
hibits a different sequence of nuclei, because 66 Ni is not 
found as an equilibrium nucleus by using the data from 
the droplet model from Myers and Swiatecki ^(| . But as 
one can see, this little deviation in the sequence of nuclei 
does not have any noticeable consequences on the equa- 
tion of state by comparing to all other graphs. A closer 
look reveals tiny jumps in the mass density for constant 
pressure. This is not surprising because such a behaviour 
is predicted by BPS [i| and explained in Sec.|Bj see Eqs. 
(|14H . Only small differences can be seen in the high den- 
sity range where the graphs begin to separate. In order 
to really see these tiny differences, we zoomed into the 
high mass density region. The jumps and the separa- 
tion of the graphs are clearly recognizable in Fig. [3 At 
high mass density, the graphs separate from each other 
because we use different nuclear models which have dif- 
ferent binding energies per nucleon for the same nuclei. 
This leads to different equilibrium nuclei and different 
equations of state. The equilibrium nuclei of known nu- 
clei are marked along the graph. The iron nucleus 56 Fe 
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is the energetically favoured one until the energy density 
reaches 10 7 g/cm 3 . Then the sequence of nuclei continues 
with the nickel isotopes 62 Ni, 64 Ni, and 66 Ni and jumps 
then to the heavier nuclei 86 Kr and 84 Se. Beyond 84 Se 
and a energy density of about 10 10 g/cm 3 , we find that 
the sequence depends on the nuclear model and the nu- 
clear mass table used. 

In Fig. |31 we show the adiabatic index T as a func- 
tion of the mass density. As one can see, there are no 
noticeable differences between the graphs of the nuclear 
models used. At high mass density, T asymptotically ap- 
proaches the value of the relativistic limit, i.e. T = |. 
At the transition points from one equilibrium nucleus to 
another, the value of the adiabatic index jumps to zero. 
For simplicity, these jumps are not shown in Fig. 



B. Spherical Relativistic Models and Effects of 
Pairing 

In Fig. 0] the effect of pairing in nuclei on the neu- 
tron dripline and on the sequence of nuclei in the outer 
crust of non-accreting cold neutron stars is depicted. We 
compare relativistic nuclear models with pairing to those 
without pairing effects for spherical calculations. The 
form of the neutron dripline shows some distinct features 
(but note for the following discussion that the mass ta- 
bles for NL-Z2 and PCF1 start only at Z = 26). In par- 
ticular, there are 'plateaus' and 'walls' visible along the 
neutron dripline indicating particular strong shell effects 
for neutrons as well as protons. In addition, the neu- 
tron dripline does not continuously increase with increas- 
ing neutron number, but exhibits 'peninsulas' of particle 
stable isotopes, especially pronounced for the set NL-Z2 
(with pairing) around N — 90—100 and Z — 32. Hence, a 
sequence of stable isotopes with increasing neutron num- 
ber stops first at the neutron dripline followed by a region 
of unstable isotopes, but then exists again a sequence of 
stable isotopes for higher N. These features were only 
found because the calculations did not stop at the neu- 
tron dripline but were performed even for several nuclei 
behind the neutron dripline. We also note, that the effect 
of pairing on the neutron dripline is rather small. The 
inclusion of pairing does not shift the neutron dripline 
but smoothes out the neutron-driplinc. 

The effect of pairing on the sequence of nuclei is also 
rather small. Differences in the equilibrium nuclei can 
be only seen at large mass numbers A. For low mass 
numbers, there is no difference because we use consis- 
tently the data of the atomic mass table [25| whenever 
available. However, pairing has an effect on the sequence 
of neutron numbers as it splits those of the point cou- 
pling model PCF1 and the chiral model: from 82 to 80 
and 82 in the point coupling model PCF1 and from 70 
to 68 and 70 in the chiral model. The reason is that 
the effect of pairing leads to the occupation of extra en- 
ergy levels which are preferred in comparison to the ones 
without pairing. The smearing of energy levels also cause 



the smoothing effect on the neutron dripline. Note, that 
the chiral model does not have the usual magic neutron 
number 82 but the magic neutron number 70. This ef- 
fect might be related to the inclusion of tensor terms in 
the chiral model which are absent in the other relativistic 
models used here. 



C. Non-relativistic Parameterizations 

In Figs. and we compare non-relativistic Skyrme- 
type models. There are no big differences visible for the 
neutron dripline and the sequence of nuclei between the 
parameter sets SkSC4, which uses the ETFSI method, 
MSk7, a Hartree-Fock calculation with BCS pairing, and 
BSk8, the full HFB calculation, in Fig. The neutron 
driplines of all models in Fig. are next to each other 
and the sequence of nuclei are almost the same. Only 
the set BSk8 has a shift to larger proton numbers Z in 
the sequence of nuclei, AZ = 2. All parameterizations 
shown in Fig. have magic neutron numbers 50 and 82. 
On the other hand, the nuclear models of Dobaczewski et 
al. show pronounced differences in the neutron dripline 
as well as in the sequences of nuclei, see Fig. The neu- 
tron dripline of the set SkP is shifted to smaller proton 
numbers compared to the set SLy4. The neutron dripline 
of the set SkM* is shifted to even smaller proton numbers 
in comparison to the neutron dripline of the set SkP. The 
sequence of nuclei of the set SLy4 is shifted to larger pro- 
ton numbers in comparison to the other two sets, with 
AZmax = 4. The sequence of nuclei for SkP makes an un- 
usual big jump to the nucleus 86 Fe and has also unusual 
magic proton numbers 30 and 38 while the other two 
sets only have the magic proton number 28. However, 
all models in Fig. have magic neutron numbers 50 and 
82. We conclude, that the details of the approximation, 
ETFSI, BCS, or HFB, are not important for the loca- 
tion of dripline and the sequence of neutron-rich nuclei 
in neutron star matter, as long as the parameters are fit- 
ted to a similar (extended) set of observables. There are, 
however, substantial differences if different sets of observ- 
ables are used for the fitting procedure. Note, that all 
mass tables used here include effects from deformations. 



D. Neutron Driplines 

We now compare relativistic and non-relativistic pa- 
rameter sets and mass tables in detail, delineating in 
particular the role of deformations. In Fig. [7| we show 
the neutron driplines of all parameter sets as listed in 
Table J] including the relativistic parameterizations NL3 
and TMA within a calculation including deformations. 

By comparing the relativistic models with and with- 
out pairing, one again recognizes that pairing makes the 
neutron driplines smoother. The prominent proton num- 
ber found in the dripline of the point coupling model 
PCF1 with pairing is 40. Its neutron numbers are 68, 
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80, and 110 with pairing, and 70, 82, and 112 without 
pairing. The corresponding proton numbers of the chiral 
model are 24, 42, and 46 with pairing, and 24, 38, 42, 
and 46 without pairing. Its neutron numbers are 60, 68, 
and 110 with pairing, and 62, 70 and 112 without pair- 
ing. The proton number found for the model NL-Z2 with 
and without pairing is 32. Its neutron numbers are 68, 
82, 102, and 112 with pairing, and 70, 82, 102, and 112 
without pairing. NL3 with pairing has many distinctive 
steps in the dripline, the proton numbers are 26, 32, 38 
and 40. The neutron numbers also for NL3 without pair- 
ing are 62, 70, 82 and 112. By comparing spherical with 
deformed nuclei, one finds that the effect of deforma- 
tion is that the neutron dripline rises steeper and nearly 
linear. The neutron driplines of the relativistic nuclear 
models with deformations, NL3 and TMA, extend from 
(Z, N) ~ (20, 45) to (Z, N) ~ (50, 105) in a nearly linear 
and direct way in stark contrast to the wiggly neutron 
driplines of the spherical relativistic calculations. The set 
TMA has less pronounced neutron and proton numbers 
in the neutron dripline in comparison to the spherical 
calculations. Noticeable occupied proton numbers of the 
model TMA are 22, 26, and 40 which is unusual, and 
a magic neutron number is 82. In some cases, the neu- 
tron dripline steps down before continuing upwards for 
the mass table of set TMA (at Z = 30, N — 74 and at 
Z = 38, N = 88). The model NL3 with deformations has 
the proton numbers 30 and 44, the neutron numbers are 
54, 82 and 86. Compared to the spherical calculations, 
the deformed calculations of NL3 and TMA show a rather 
similar straight behaviour, although the precise endpoint 
of a certain sequence of isotopes can differ drastically, in 
particular for Z = 28. It is evident, that the inclusion of 
deformations is crucial for the overall shape of the neu- 
tron dripline but that the location depends strongly on 
the parameterization. 

We now discuss the (deformed) Skyrme calculations in 
comparison. All neutron driplines of the nuclear mod- 
els taken from BRUSLIB do not show big differences. 
Only the neutron dripline of the model SkSC18 decreases 
sometimes to much lower proton numbers Z and then 
quickly rises again. We attribute this effect to the ap- 
proximate treatment of shell corrections to the bind- 
ing energy which will be particular important close to 
the dripline. All neutron driplines extracted from the 
BRUSLIB mass tables range from (Z, N) ~ (20, 48) to 
(Z, N) ~ (50, 108) in a nearly linear fashion. A notice- 
able proton number in the models of BRUSLIB is 30 and 
the neutron number 82. By comparing the nuclear mod- 
els from BRUSLIB with the ones from Dobaczewski et al., 
one can see that there are no substantial differences be- 
tween the location of the neutron driplines. Only the neu- 
tron driplines of the models SkM* and SkP are shifted to 
lower proton numbers compared to the other parameter 
sets. The neutron driplines of the sets SLy4 and SLy4HO 
(set SLy4 calculated in hydrogen oscillator basis) ex- 
tends from (Z, N) ~ (20, 48) to (Z, N) ~ (50, 108) while 
that one of the set SkM* goes from (Z, N) = (20, 56) to 



(Z, N) = (46, 115), and the one from the model SkP from 
{Z,N) = (20,46) to (Z,N) = (50,114). By comparing 
the neutron driplines of the models from Dobaczewski et 
al. to those from BRUSLIB, one recognizes that there 
are much more marked shells effects in the models from 
Dobaczewski et al. Noticeable proton numbers of the 
models SLy4 and SLy4HO are 32, 42, and 46 and the 
neutron number 82. A marked proton shell of the set 
SkM* is 30. A marked neutron shell of this set is 82. 
The set SkP exhibits the unusually proton number 22 in 
its dripline, like the relativistic set TMA. The set SkM* 
exhibits a noticeable downward sequence of the neutron 
dripline around Z = 40 and N = 100. Similar features of 
downward shifts in the neutron dripline are seen in the 
sequences of the parameter sets from BRUSLIB around 
N = 104 and N = 108. We checked that these features 
are also present for the HFB set BSk8 and not due to the 
approximation scheme used (BCS or ETFSI). 

The neutron dripline of the finite range droplet model 
FRDM ranges from (Z,N) = (20,50) to (Z,N) = 
(50,112). It has four pronounced proton numbers: 28, 
34, 44 and 48. A noticeable neutron shell is again 82. 
By comparing all neutron driplines from Fig. {7\ one rec- 
ognizes that all non-relativistic models show the magic 
neutron number 82 to be present at the neutron dripline 
and that all proton shell closures of the non-relativistic 
models are less noticeable in comparison to it. These 
findings are different to the ones for the relativistic mod- 
els (even when including deformations) in Fig. [7| where 
both, proton and neutron shells, are more noticeable. 



E. Sequences of Nuclei 

In Fig- El we show the sequences of nuclei in the outer 
crust of non-accreting cold neutron stars by using various 
nuclear mass tables. By comparing the sequences of nu- 
clei with pairing to those without pairing, one recognizes 
that pairing has the effect of splitting the magic neutron 
numbers. The magic neutron number 82 is split into 
80 and 82 in the chiral model and in the point coupling 
model PCF1. By comparing spherical with deformed cal- 
culations, one can see that the effects of deformation on 
the sequence of nuclei is small. Models with and without 
deformation have the same magic neutron numbers: 50 
and 82. They also have the same sequence of neutron- 
rich nuclei in the nuclear chart. An exceptional case is 
the sequence of nuclei of the chiral model which shows 
differences concerning the magic neutron numbers. The 
sequence of nuclei within the chiral model does not fol- 
low the magic neutron number 82 but instead continues 
down the neutron number 70 basically all the way to 
the neutron dripline. Again, we attribute this behaviour 
to the inclusion of tensor terms for the p meson in the 
effective Lagrangian which induce a different isospin de- 
pendent spin-orbit strength compared to the other rela- 
tivistic models. 

The different sequences of nuclei for the parameter sets 
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of BRUSLIB are nearly on top of each other and are 
strikingly along the magic neutron numbers 50 and 82. 
There a only slight differences, as the sequence for set 
BSk8 reaches Z = 46 at N = 82, while all other sets 
from BRUSLIB start the sequence along N = 82 at Z = 
44. The set SLy4 features similar strong correlations in 
the sequence of nuclei along TV = 50, 82, but like the 
relativistic sets NL3 and TMA populates at maximum 
Z = 42 for N = 82. 

By comparing the sets from BRUSLIB and SLy4 to 
those of SkM* and SkP one recognizes that the sequences 
of nuclei from SkM* and SkP are distinctly different in 
their paths in the nuclear chart as they mainly follow 
along isotopes and not along the magic neutron numbers 
50 and 82. This finding is in stark contrast to the other 
parameterizations and mass tables used here, making the 
sets SkP and SkM* noticeable exceptions from our gen- 
eral results. 

In particular, the more modern Skyrme HFB mass ta- 
bles of the sets SLy4 and BSk8 do not exhibit a pro- 
nounced sequence along isotopes but are morestuck to 
the magic neutron numbers 50 and 82. The general trend 
of the sequence of nuclei being along these magic neutron 
numbers is supported by the calculation for the FRDM 
which is strikingly similar to the results for the sets TMA 
and SLy4. Surprisingly, the classic sequence of BPS, al- 
though using rather old mass tables, arrived at a similar 
sequence for large mass numbers compared to the most 
modern mass tables! 



the sequences of nuclei in /3-equilibrium of the modern 
sets shown do not exhibit pronounced differences, on the 
contrary, the sequences are nearly on top of each other. 
There are only small exceptions, as there is a shift to 
larger proton numbers in the sequence of nuclei for the 
set BSk8 with AZ max = 4. The sets BSk8 and NL3 have 
one nucleus in their sequence with N = 52, set NL3 ends 
its sequence at N = 84. Besides that and the starting 
iron and nickel isotopes, the sequences of nuclei for all the 
selected models and mass tables follow tenaciously the 
magic neutron numbers 50 and 82 throughout the (un- 
known) nuclear chart until hitting the neutron dripline. 

The endpoint of the sequence of nuclei for all the mod- 
ern mass tables studied here is quite similar and hap- 
pens to be at Z = 34 - 38 with N = 82 (N = 84 for 
the set NL3). In this region of the nuclear chart, the 
neutron driplines of the various sets examined here do 
also demonstrate to be rather similar. These features 
most likely originate from the fact that the magic neu- 
tron number 82 prevails up to the neutron dripline in the 
nuclear mass tables used here, which however, relies on a 
substantial extrapolation from the masses of known nu- 
clei up to the neutron dripline. Shell quenching effects, as 
advocated by HZD |17| . can possibly exist, in particular 
for different isospin dependences of the spin-orbit terms. 
Then the sequence of nuclei as well as its endpoint can be 
located quite differently as seen e.g. for sets SkM*, SkP 
and Chiral. 



F. Neutron Driplines and Sequences of Nuclei of 
State-of-the-Art Nuclear Models 

In Fig. El we compare the most modern and state-of- 
the-art mass tables of all the sets listed in Tabled There 
are no drastic differences in the shape of the neutron 
dripline. The neutron driplines of all the sets presented 
in the figure have an approximately linear behaviour ex- 
tending from (Z,N) ~ (20,46) to (Z, N) = (36,82). At 
the magic neutron number 82, the neutron driplines sud- 
denly change their slope and follow a vertical path from 
(Z,N) = (36,82) to (Z,N) = (40,82). Note that at the 
magic neutron number 82, the neutron driplines of all the 
selected sets are equal. Therefore, the region at N = 82 
and around Z = 38 shows to be the one with the smallest 
difference in the precise location of the neutron dripline 
for all modern sets used here. From (Z, N) = (40, 82) 
to (Z.N) ~ (50,108), the neutron driplines again have 
approximately the same gradient as for low mass num- 
bers. But there are differences in the regions of the nu- 
clide chart where Z ~ 28 and N ~ 62: If one compares 
the neutron driplines from the mass tables of sets TMA 
and NL3 (which includes effects from deformations), then 
one can see that their neutron driplines differ by AZ = 4 
and AN = 16. There are also differences of similar range 
for the neutron driplines at (Z, N) ~ (44, 100) if one 
compares the model TMA with SLy4 and NL3. Also, 



Finally, we note that in all our calculations, there do 
not appear any super-heavy nuclei in the sequences of 
nuclei of the nuclear models and sets which we use in 
this work. We checked that superheavy elements can 
indeed appear in the sequence of neutron-rich nuclei in 
the outer crust of neutron stars if we artificially increase 
their binding energies per nucleon b by about one MeV. 

In Fig. 1101 we show the dependence of the proton num- 
ber Z on the mass density p by using the representative 
models of Table [I] In this plot, we restrict ourselves to 
the high mass density region p > 10 10 g/cm 3 because for 
10 4 g/cm 3 < p < 3.5 x 10 10 g/cm 3 we obtain the same 
elements for each set as they are fixed by the experimen- 
tal nuclear data of [2i|: 26Fe, 2sNi, 36Kr, 34Se, 32 Ge, and 
30 Zn. For larger mass densities, the proton numbers of 
most of the sets shown are the same: 2sNi, 42M0, 4oZr, 
38Sr, and 36Kr. The only differences for these mass tables 
is that they arrive at the corresponding proton numbers 
at different mass densities. For the set TMA, 34Se and 
not 36Kr appears at the highest mass densities in the 
outer crust of non-accreting cold neutron stars. The set 
BSk8 is the only one that shows significant differences in 
comparison to the other sets plotted: For mass densities 
of p ~ 10 11 g/cm 3 , the BSk8 model has quite different 
proton numbers as 46Pd and 44RU appear in the sequence 
of nuclei. 
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G. Comparison to previous works 

In order to compare the sequences of nuclei of BPS 0, 
HZD i and HP 20] with our new results as presented 
here, we list all the sequences of nuclei of these previous 
works in tabular form. Table El shows the sequence of 
nuclei of BPS, Tables EH and EH the sequences of nu- 
clei of HZ D 1171 , Tables and E] the sequences of nu- 
clei of HP [2~of7^ Finally. Tables IVTT1 and lyTTTl summarize 
the sequences of nuclei obtained in this work by using 
the modern nuclear mass tables from BSk8 and TMA as 
the characteristic ones of our whole set of mass tables 
investigated. Note that BPS and HZD only use theoret- 
ically computed mass tables, while we, as HP, incorpo- 
rate experimental data tables in addition. As we know 
from Fig. the equation of state is not affected signifi- 
cantly by a different sequence of nuclei. That is why we 
do not show the number density, the pressure and the 
mass density in every table. Besides that, we are mainly 
interested in the sequence of nuclei here. As one can 
see from the tables, the nucleus 66 Ni is not present in 
the crust of non-accreting cold neutron stars if one uses 
the droplet model from Myers and Swiatecki 0] or the 
spherical model SkM* from Dobaczewski, Flocard, and 
Treiner [3. By taking newer experimental data, this 
nucleus always exists in the crust. Even with the droplet 
model of Myers 0] , this nucleus is found. A second dif- 
ference in the sequence of nuclei of BPS compared to 
HZD, HP, BSk8, and TMA is that in the latter three one, 
76 Fe never occurs. It only occurs by using the droplet 
model from Myers and Swiatecki or that one from 
Myers When we compare the sequence of nuclei of 
the model BSk8 with the sequence of nuclei from HZD, 
HP and that from BPS, we recognize that the nucleus 
78 Ni is not present in the sequence of nuclei of the model 
BSk8. The reason why 78 Ni appears in the sequence of 
nuclei from HZD and not in the sequence of nuclei of the 
BSk8 model is that HZD used the atomic mass table from 
1993 [2ll and we use the newer one from 2003 Hj| . With 
the set TMA, one obtains a similar sequence of nuclei as 
for the droplet model of Moller and Nix (2^ , see TablesM 
and IVIlJ The only differences are that by using the set 
TMA, 126 Ru is not present in the sequence of nuclei and 
that the sequence ends with the nucleus 116 Se and not 
with 118 Kr. 

The sequence of nuclei obtained by HZD 01 using the 
set SkP within a spherical calculation of Dobaczewski, 
Flocard, and Treiner (see table llllfl is drastically dif- 
ferent from all other sequences of nuclei listed here. Un- 
usual nuclei for the crust, such as 68 Ni or later in the se- 
quence, nickel nuclei with large neutron numbers appear. 
In addition, our results for the same set SkP differ from 
the results of HZD 17] . Note, that the mass table used 
by HZD is based on a spherical HFB calculation (and a 
different pairing force) while the one used in this work 
includes effects from deformations and a modern (5-force 
pairing. The sequences of nuclei for both mass tables fol- 
low closely a certain proton number for a wide range of 



isotopes, but HZD find Z = 28 while we find Z = 30 to 
be prominent. Only the set SkM* shows a similarly wide 
sequence of nuclei with constant Z, however appearing 
for Z = 28 and not for Z = 30. The sequences of nu- 
clei found here differ from the one by HZD also at large 
mass numbers. HZD find, that the sequence jumps from 
Z = 28 to Z = 40 ending with the nucleus 134 Zr while 
the sequence in our calculation populates Z — 38 (and 
once Z = 26) and ends with the nucleus 126 Sr. We stress, 
that we compare here different mass tables which are, al- 
though based on the same parameter set SkP, computed 
in significantly different approximations (spherical versus 
deformed) and use different pairing forces. The general 
trend of the sequence of neutron-rich nuclei, however, is 
rather similar, as the sequence follows mainly isotope se- 
quences (constant Z) and not the magic neutron numbers 
50 and 82. 

The last nucleus in the sequences of nuclei found by 
using different theoretical nuclear models are shown in 
Table |IX| All of the models listed in Table ED except 
NL3 and that of 0] have the magic neutron number 82. 
The charge of the final nucleus of the sequence varies 
between Z = 32 — 40 for previous works while we find a 
range of Z — 34 — 38 for the modern mass tables. 



V. SUMMARY AND CONCLUSIONS 

In this work, we investigated the outer crust of non- 
accreting cold neutron stars by using the BPS model and 
the newest nuclear data. If data from the atomic mass 
table of 2003 from Audi, Wapstra, and Thibault [2^] was 
available for the corresponding nuclei, we always pre- 
ferred using this experimental data. We updated pre- 
vious work of BPS H, HZD [HI, and HP HJ which were 
based on older data and/or mass parameterizations. We 
also compared various nuclear models as listed in Ta- 
ble H] in order to check their differences concerning the 
neutron dripline, magic neutron numbers, the equation 
of state, and the sequence of nuclei in the outer crust 
of non-accreting cold neutron stars. To our knowledge, 
this is the first detailed investigation of this kind for an 
enlarged set of most modern nuclear models and state- 
of-the-art theoretical mass tables, which entails pairing 
effects, includes effects of deformation and studies rela- 
tivistic models in comparison to non-relativistic ones. 

We obtain the following results: The equation of state 
is not affected by small differences in the sequence of nu- 
clei. There are jumps in the mass density at constant 
pressure if the equilibrium nucleus of the ground state 
changes to another one. The adiabatic index T is not 
affected by changes in the sequence of nuclei and by us- 
ing data from different nuclear models. The location of 
the neutron driplines of the Skyrme type models SkP 
and SkM* are at smaller proton numbers Z compared to 
the other Skyrme type models used here (SLy4 and the 
ones of BRUSLIB). The neutron dripline rises steeper 
and is nearly linear in mass tables including effects from 
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deformation as compared to spherical calculations. Pair- 
ing has the effect of smoothing out steps on the neu- 
tron dripline found when pairing effects are switched off. 
Besides that, pairing has an impact on the magic neu- 
tron numbers towards the neutron driplines. It splits the 
magic neutron number of the point coupling model PCF1 
from 82 to 80 and 82 and the magic neutron number of 
the chiral model from 70 to 68 and 70. The sequences 
of nuclei obtained by using theoretical nuclear models 
can have some differences at large mass numbers A be- 
cause they predict different binding energies. For low 
mass numbers, as we use the experimental data from the 
atomic mass table, the sequences of nuclei in the outer 
crust of neutron stars are found to be identical. The 
nucleus 66 Ni does not occur in the sequences of nuclei 
of BPS and HZD but it is present in the sequences of 
nuclei if one uses newer nuclear data. The sequence of 
nuclei usually found in all parameterizations used pro- 
ceeds along the magic neutron numbers 50 and 82 all 
the way to the neutron dripline. Only the sets for SkP 
and SkM* show sequences of nuclei which follow a cer- 
tain (magic) proton number, i.e. Z = 28, 30 and Z = 38, 
which is in contrast to all other sets investigated here, 
being non-relativistic or relativistic ones. 

By comparing a selected set of modern models and 
mass tables, the location of the neutron dripline as well 
as the sequence of nuclei in /3-equilibrium are found to be 
strikingly similar. The sequence of nuclei follows nearly 
entirely the magic neutron numbers 50 and 82. The set 
BSk8 differs slightly from all the other modern models as 
the sequence of nuclei is shifted to larger proton numbers 
while following the standard path along N — 50 and 
82. The endpoint of the sequence coincides with a region 
in the nuclear chart where the location of the neutron 
dripline of the various models investigated here coincides. 
Hence, the final nucleus in the sequence can be pinned 
down in a rather narrow range and is extracted to be 
around Z = 36 and N = 82. 

Finally, we note that in all our calculations, there do 
not appear any super-heavy nuclei in the sequences of 
nuclei of the nuclear models which we use in this paper. 
Some of them occur if the binding energies per nucleon b 



is increased artificially by about one MeV. 

The results presented here rely only on the binding 
energy or mass of neutron-rich isotopes in the range of 
Z = 26 to Z = 40, a region of the nuclear chart which 
will be explored in the near future by the FAIR facility at 
GSI, Darmstadt, as well as by the ISAC facility at TRI- 
UMF and by the RIA project. With these measurements, 
the non-relativistic and relativistic models utilized in this 
work will be tightly constrained on their isospin proper- 
ties towards the neutron dripline and one can shed more 
light on the actual sequence of neutron-rich nuclei in the 
outer crust of neutron stars. We point out that the pro- 
ton number of nuclei is important for heat transport and 
for Ohmic dissipation of the magnetic field. Other ex- 
traterrestrial sources of information might come from the 
observations of x-ray binary systems, where a neutron 
star is accreting material from its stellar companion. X- 
ray bursts are generated by accretion in the outer crust 
of neutron stars which transport material of the crust to 
the surface. The x-ray satellites XMM-Newton, Chandra 
or in the near future XEUS can give a detailed spectro- 
scopic a naly sis of an x-ray burst. Cottam, Paerels, and 
Mendez |4|j reported about redshifted spectral lines of 
highly ionized elements in x-ray bursts of the low-mass 
x-ray binary EXO 0748-676. 
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FIG. 1: The equations of state, i.e. the pressures as a function of the mass densities, calculated by using the BPS model and 
the binding energies of various nuclear models (NL3p: set NL3, spherical calculation with pairing, PCFlnp: set PCF1 without 
pairing, NL3def: set NL3 in the deformed calculation). 



H [MeV] 


A» e [MeV] 


Pmax [g/cm 3 ] 


P [dy 


nc 


/cm 2 ] 


rib 


cm 3 ] 


Element 


Z 


N 


930.60 


0.95 


8.09 x 10 6 


5.29 


X 


10 23 


4 


88 


X 


10 30 


56 Fe 


26 


30 


931.31 


2.60 


2.69 x 10 8 


6.91 


X 


10 25 


1 


62 


X 


10 32 


62 Ni 


28 


34 


932.00 


4.24 


1.24 x 10 9 


5.20 


X 


10 26 


7 


48 


X 


10 32 


64 Ni 


28 


36 


933.33 


7.69 


8.15 x 10 9 


5.78 


X 


10 27 


4 


90 


X 


10 33 


84 Se 


34 


50 


934.42 


10.61 


2.23 x 10 10 


2.12 


X 


10 28 


1 


34 


X 


10 34 


82 Ge 


32 


50 


935.48 


13.58 


4.88 x 10 10 


5.70 


X 


10 28 


2 


93 


X 


10 34 


80 Zn 


30 


50 


937.68 


19.97 


1.63 x 10 11 


2.68 


X 


10 29 


9 


74 


X 


10 34 


78 Ni 


28 


50 


937.78 


20.25 


1.78 x 10 11 


2.84 


X 


10 29 


1 


07 


X 


10 35 


76 Fe 


26 


50 


937.83 


20.50 


1.86 x 10 11 


2.93 


X 


10 29 


1 


12 


X 


10 35 


124 Mo 


42 


82 


938.57 


22.86 


2.67 x 10 11 


4.55 


X 


10 29 


1 


60 


X 


10 35 


122 Zr 


40 


82 


939.29 


25.25 


3.73 x 10 11 


6.79 


X 


10 29 


2 


23 


X 


10 35 


120 gr 


38 


82 


939.57 


26.19 


4.32 x 10 11 


7.87 


X 


10 29 


2 


59 


X 


10 35 


118 Kr 


36 


82 



TABLE II: Sequence of nuclei in the outer crust of non-accreting cold neutron stars calculated by using the nuclear data of 
Myers and Swiatecki fl^ to reproduce the results of BPS The last line corresponds to the neutron drip point. 
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FIG. 2: The same as Fig. □ for the high mass density range. 



Me [MeV] 


pmax 


[g/cm 3 ] 


Element 


Z 


N 


0.95 


7.96 x 10 6 


56 Fe 


26 


56 


2.61 


2.70 x 10 8 


62 Ni 


28 


34 


4.17 


1.18 x 10 9 


64 Ni 


28 


36 


6.94 


5.88 


X 


10 9 


68 Ni 


28 


40 


9.12 


1.36 


X 


10 10 


84 Se 


34 


50 


9.16 


1.40 


X 


10 10 


70 Ni 


28 


42 


10.06 


1.93 


X 


10 10 


72 Ni 


28 


44 


10.86 


2.43 


X 


10 10 


78 Zn 


30 


48 


13.24 


4.52 


X 


10 10 


80 Zn 


30 


50 


13.58 


4.97 


X 


10 10 


76 Ni 


28 


48 


16.51 


9.17 


X 


10 10 


78 Ni 


28 


50 


18.04 


1.23 


X 


10 11 


80 Ni 


28 


52 


19.38 


1.56 


X 


10 11 


82 Ni 


28 


54 


21.12 


2.07 


X 


10 11 


84 Ni 


28 


56 


22.89 


2.70 


X 


10 11 


86 Ni 


28 


58 


24.36 


3.33 


X 


10 11 


88 Ni 


28 


60 


25.10 


3.77 


X 


10 11 


130 Zr 


40 


90 


25.97 


4.25 


X 


10 11 


132 Zr 


40 


92 


26.26 


4.45 


X 


10 11 


134 Zr 


40 


94 



TABLE III: Sequence of nuclei in the outer crust of non-accreting cold neutron stars as calculated by Haensel, Zdunik, and 
Dobaczewski Jrjfl ■ Nuclear masses are taken from a spherical calculation using the parameter set SkP of Dobaczewski, Flocard, 
and Treiner 11811. 
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FIG. 3: The adiabatic index as a function of the mass density, calculated by using the BPS model and the binding energies of 
various nuclear models. The value of Y of the horizontal line is equal to = . 



Me [MeV] 


Pmax [g/cm 3 ] 


Element 


Z 


N 


0.95 


7.96 x 10 6 


56 Fe 


26 


30 


2.61 


2.70 x 10 8 


62 Ni 


28 


34 


4.28 


1.29 x 10 9 


64 Ni 


28 


36 


4.57 


1.61 x 10 9 


66 Ni 


28 


38 


5.32 


2.63 x 10 9 


68 Ni 


28 


40 


6.21 


4.34 x 10 9 


80 Ge 


32 


48 


9.69 


1.70 x 10 10 


82 Ge 


32 


50 


12.26 


3.59 x 10 10 


80 Zn 


30 


50 


18.22 


1.23 x 10 11 


78 Ni 


28 


50 


18.73 


1.41 x 10 11 


76 Fe 


26 


50 


20.15 


1.83 x 10 11 


122 Zr 


40 


82 


22.19 


2.53 x 10 11 


i20 Sr 


38 


82 


24.24 


3.42 x 10 11 


118 Kr 


36 


82 


26.28 


4.55 x 10 11 


116 Se 


34 


82 


26.82 


5.05 x 10 11 


114 Go 


32 


82 



TABLE IV: Sequence of nuclei in the outer crust of non-accreting cold neutron stars calculated by Haensel, Zdunik, and 
Dobaczewski flTf . Nuclear masses are taken from the droplet model of Myers [l9|. 
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FIG. 4: Nuclear charts of various relativistic nuclear models and their neutron driplines. The upper plots correspond to nuclear 
models with pairing while the lower ones correspond to nuclear models without pairing. The crosses mark the nuclei which 
are taken from the atomic mass table [2^|. The thick lines and the points show the sequence of nuclei in the outer crust of 
non-accreting cold neutron stars by using various nuclear data. 




FIG. 5: Nuclear chart of various non-relativistic deformed Skyrme Hartree-Fock-Bogolyubov calculations as taken from BRUS- 
LIB |2(J [53 and their neutron driplines. The crosses mark the nuclei which are taken from the atomic mass table (25ft . The 
thick lines and points show the sequence of nuclei in the outer crust of non-accreting cold neutron stars by using various nuclear 
models. 
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FIG. 6: The same as Fig. but now for non-relativistic deformed Skyrme Hartree-Fock-Bogolyubov models taken from 
Dobaczewski et al. [2§|| . 




FIG. 7: Nuclear charts showing the neutron driplines of various theoretical nuclear models (see TableQJ. The crosses mark the 
nuclei which are taken from the atomic mass table Hal. 
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FIG. 8: Same as Fig. [7|but plotting here the sequence of nuclei in the outer crust of non-accreting cold neutron stars by using 
various nuclear models. 




FIG. 9: Nuclear chart of the selected most modern models from Table Q and their neutron driplines. The crosses mark the 
nuclei which are taken from the atomic mass table |25l| . The thick lines and points show the sequence of nuclei in the outer 
crust of non-accreting cold neutron stars. 
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FIG. 10: The proton number Z as a function of the mass density p for the selected most modern models as depicted also in 
Fig. 



fie [MeV] 


Pmax [g/cm 3 ] 


Element 


Z 


N 


0.95 


7.96 x 10 6 


56 Fe 


26 


30 


2.61 


2.71 x 10 8 


62 Ni 


28 


34 


4.31 


1.30 x 10 9 


64 Ni 


28 


36 


4.45 


1.48 x 10 9 


66 Ni 


28 


38 


5.66 


3.12 x 10 9 


86 Kr 


36 


50 


8.49 


1.10 x 10 10 


84 Se 


34 


50 


11.44 


2.80 x 10 10 


82 Ge 


32 


50 


14.08 


5.44 x 10 10 


80 Zn 


30 


50 


16.78 


9.64 x 10 10 


78 Ni 


28 


50 


18.34 


1.29 x 10 11 


126 Ru 


44 82 


20.56 


1.88 x 10 11 


124 Mo 


42 


82 


22.86 


2.67 x 10 11 


122 Zr 


40 


82 


25.38 


3.79 x 10 11 


120 Sr 


38 


82 


26.19 


4.33 x 10 11 


118 Kr 


36 


82 



TABLE V: Sequence of nuclei in the outer crust of non-accreting cold neutron stars calculated by Haensel and Pichon [2(J. 
Upper part: using experimental nuclear masses. Lower part: using binding energies from the mass formula of Moller and 
Nix |22J. The last line corresponds to the neutron drip point. 
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Me [MeV] 


Pmax [g/cm 3 ] 


Element 


Z 


N 


17.44 


1.08 x 10 11 


78 Ni 


28 


50 


19.13 


1.47 x 10 11 


126 Ru 


44 82 


21.66 


2.20 x 10 11 


124 Mo 


42 


82 


24.13 


3.15 x 10 11 


122 Zr 


40 


82 


26.05 


4.10 x 10 11 


120 Sr 


38 


82 



TABLE VI: Sequence of nuclei in the outer crust of non-accreting cold neutron stars calculated by Haensel and Pichon [2fT |. 
Upper part: last experimental nuclear mass. Lower part: using the mass formula of Aboussir et al. ,231. The last line 
corresponds to the neutron drip point. 



[i [MeV] 


/i e [MeV] 


Pmax [g/Cm ] 


P [dyne/ cm 


r — 8i 
n b [cm J 


Element 


Z 


N 


930.60 


0.95 


8.02 x 10 6 


5.22 x 10 23 


4.83 x 10 30 


56 Fe 


26 


30 


931.32 


2.61 


2.71 x 10 8 


6.98 x 10 25 


1.63 x 10 32 


62 Ni 


28 


34 


932.04 


4.34 


1.33 x 10 9 


5.72 x 10 26 


8.03 x 10 32 


64 Ni 


28 


36 


932.09 


4.46 


1.50 x 10 9 


6.44 x 10 26 


9.04 x 10 32 


66 Ni 


28 


38 


932.56 


5.64 


3.09 x 10 9 


1.65 x 10 27 


1.86 x 10 33 


86 Kr 


36 


50 


933.62 


8.38 


1.06 x 10 10 


8.19 x 10 27 


6.37 x 10 33 


84 Se 


34 


50 


934.75 


11.43 


2.79 x 10 10 


2.85 x 10 28 


1.68 x 10 34 


82 Ge 


32 


50 


935.89 


14.61 


6.07 x 10 10 


7.63 x 10 28 


3.65 x 10 34 


80 Zn 


30 


50 


936.44 


16.17 


8.46 x 10 10 


1.15 x 10 29 


5.08 x 10 34 


82 Zn 


30 


52 


936.63 


16.81 


9.67 x 10 10 


1.32 x 10 29 


5.80 x 10 34 


128p d 


46 


82 


937.41 


19.16 


1.47 x 10 11 


2.23 x 10 29 


8.84 x 10 34 


126 Ru 


44 


82 


938.12 


21.35 


2.11 x 10 11 


3.45 x 10 29 


1.26 x 10 35 


124 Mo 


42 


82 


938.78 


23.47 


2.89 x 10 11 


5.05 x 10 29 


1.73 x 10 35 


122 Zr 


40 


82 


939.47 


25.77 


3.97 x 10 11 


7.36 x 10 29 


2.38 x 10 35 


120 gr 


38 


82 


939.57 


26.09 


4.27 x 10 11 


7.74 x 10 29 


2.56 x 10 35 


118 Kr 


36 


82 



TABLE VII: Sequence of nuclei in the outer crust of non-accreting cold neutron stars calculated by using the experimental 
nuclear data from the atomic mass table |25| (upper part), and the theoretical mass table of the Skyrme model BSk8 as listed 
by BRUSLIB (lower part). Note that the experimental data from the atomic mass table |25| is always taken first if available. 
The last line corresponds to the neutron drip point. 



H [MeV] 


A» e [MeV] 


Pmax [g/cm 3 ] 


P [dyne/cm 2 ] 


n b [cm 3 ] 


Element 


Z 


N 


930.60 


0.95 


8.02 x 10 6 


5.22 x 10 23 


4.83 x 10 30 


56 Fe 


26 


30 


931.32 


2.61 


2.71 x 10 8 


6.98 x 10 25 


1.63 x 10 32 


62 Ni 


28 


34 


932.04 


4.34 


1.33 x 10 9 


5.72 x 10 26 


8.03 x 10 32 


64 Ni 


28 


36 


932.09 


4.46 


1.50 x 10 9 


6.44 x 10 26 


9.04 x 10 32 


66 Ni 


28 


38 


932.56 


5.64 


3.09 x 10 9 


1.65 x 10 27 


1.86 x 10 33 


86 Kr 


36 


50 


933.62 


8.38 


1.06 x 10 10 


8.19 x 10 27 


6.37 x 10 33 


84 Se 


34 


50 


934.75 


11.43 


2.79 x 10 10 


2.85 x 10 28 


1.68 x 10 34 


82 Ge 


32 


50 


935.93 


14.71 


6.21 x 10 10 


7.86 x 10 28 


3.73 x 10 34 


80 Zn 


30 


50 


937.28 


18.64 


1.32 x 10 11 


2.03 x 10 29 


7.92 x 10 34 


78 Ni 


28 


50 


937.63 


19.80 


1.68 x 10 11 


2.55 x 10 29 


1.01 x 10 35 


124 Mo 


42 


82 


938.13 


21.38 


2.18 x 10 11 


3.48 x 10 29 


1.31 x 10 35 


122 Zr 


40 


82 


938.67 


23.19 


2.89 x 10 11 


4.82 x 10 29 


1.73 x 10 35 


120 gr 


38 


82 


939.18 


24.94 


3.73 x 10 11 


6.47 x 10 29 


2.23 x 10 35 


118 Kr 


36 


82 


939.57 


26.29 


4.55 x 10 11 


8.00 x 10 29 


2.72 x 10 35 


116 Se 


34 


82 



TABLE VIII: The same as Table IVTT1 but for the relativistic nuclear model TMA. 
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Year 


Element 


Z N Model 


Refs. 


1966 


118 Kr 


36 82 Droplet model of Myers and Swiatecki 


[16] 


1977 


114 Ge 


32 82 Droplet model of Myers 


[19] 


1984 


134 Zr 


40 94 Skyrme model SkP (spherical) 


[18] 


1988 


118 Kr 


36 82 Droplet model of Moller and Nix 


[22] 


1992 


120g r 


38 82 Droplet model of Aboussir et al. 


[23] 


1995 


118 Kr 


36 82 Finite range droplet model FRDM 


[34] 


1997/1999 


120 Kr 


36 84 NL3 (with deformations) 


[44,45] 


1998/2004 


120g r 


38 82 SLy4 (with deformations) 


[28,37] 


2004 


118 Kr 


36 82 BSk8 (with deformations) 


[27] 


2005 


116 Se 


34 82 TMA (with deformations) 


[29] 



TABLE IX: The endpoint of the sequence of nuclei obtained by using different theoretical nuclear models. Upper part: results 
obtained in previous work by BPS, HZD, and HP with older nuclear models. Lower part: results obtained in this work by 
using modern models and mass tables (see Table 0. 



